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• Premise of the study: A set of microsatellite markers were developed to characterize the level of genetic diversity and gene flow 
in two plant species endemic to the Qinghai-Tibet Plateau, Spiraea alpina and S. mongolica. 

• Methods and Results: Using the Fast Isolation by AFLP of Sequences Containing repeats (FIASCO) method, 16 microsatellite 
loci showed polymorphisms in both species. In two populations of each species, the number of alleles per locus ranged from 
three to 18 in S. alpina and from four to 30 in S. mongolica. 

• Conclusions: These microsatellite markers provide an efficient tool for population genetic studies and will be used to assess the 
genetic diversity and spatial genetic structure of S. alpina and S. mongolica. 
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Spiraea alpina Pall, and S. mongolica Maxim. (Rosaceae 
subfam. Spiraeoideae) are perennial shrubs, found in western 
China and some areas of Mongolia and Siberia. The two alpine 
plants usually grow on sunny slopes or ridges. They are wide- 
spread across the Qinghai-Tibet Plateau and adjacent high- 
lands, at altitudes of 2000-4500 m (Lu et al., 2003). Due to high 
levels of morphological variation, the genus Spiraea L. has 
been classified in several ways by different authors into various 
subgenera, sections, and series (Lu et al., 2003; Potter et al., 
2007). A recent phylogeographic analysis of cpDNA variations 
in S. alpina indicated that this alpine shrub survived in multiple 
refugia during the Last Glacial Maximum and that earlier gla- 
ciations may have triggered deep intraspecific divergence 
(Zhang et al., 2012). However, the phylogeographic analysis 
based on one uniparentally inherited cpDNA fragment may 
only partly recover the phylogeographic history of a species. 
Biparentally inherited simple sequence repeat (SSR) markers 
with more polymorphism and information are necessary for 
a better understanding of the genetic structure and phylogeo- 
graphic history of S. alpina and S. mongolica. In this study, we 
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isolated 16 polymorphic microsatellite primers to facilitate the 
investigation in further studies for these two species. 

METHODS AND RESULTS 

Total genomic DNA was extracted from silica gel-dried leaves of S. alpina 
following the cetyltrimethylammonium bromide (CTAB) method (Doyle and 
Doyle, 1987). Microsatellite loci from an enriched (AG) n library were isolated 
using the Fast Isolation by AFLP of Sequences Containing Repeats (FIASCO) 
method with minor modifications (Zane et al., 2002). Approximately 300 ng 
of genomic DNA were completely digested with Msel (New England Biolabs, 
Beverly, Massachusetts, USA), and then ligated to an Msel AFLP adapter 
(S'-TACTCAGGACTCAT^VS'-GACGATGAGTCCTGAG-S') using T4 DNA 
ligase (New England Biolabs). The diluted digestion-ligation mixture (1 : 10) 
was amplified with adapter-specific primers (5'-GATGAGTCCTGAGTAAN-3'). 
For enrichment, the PCR products were denatured at 95 °C for 5 min, then hy- 
bridized with two 5'-biotinylated probes, (AC) 15 and (AG) 15 , respectively, in a 
250-jiL hybridization solution (4x saline sodium citrate [SSC], 0.1% sodium 
dodecyl sulfate [SDS], 0.5 umol/L probe) at 48°C for 2 h. Streptavidin-coated 
magnetic beads (New England Biolabs) were used to separate and capture the 
DNA fragments hybridized to the probe at room temperature for 20 min, fol- 
lowed by two washing steps: three times in TEN 100 (10 mM Tris-HCl, 1 mM 
EDTA, 100 mM NaCl [pH 7.5]) for 8 min and three times in TEN 1000 (10 mM 
Tris-HCl, 1 mM EDTA, 1 M NaCl [pH 7.5]) for 8 min. The separated single- 
stranded DNA fragments were amplified with adapter- specific primers as de- 
scribed above. The PCR products, after purification using a CASpure PCR 
Purification Kit (Sangon, Shanghai, China), were ligated into the pGEM-T 
Easy Vector (Promega Corporation, Madison, Wisconsin, USA) according to 
the manufacturer's instructions, then transformed into Escherichia coli TOP 10 
competent cells (Trans Gen Biotech, Beijing, China). Transformants were 
plated, and insert-containing clones were selected by blue-white screening with 
ampicillin, X-Gal, and isopropyl-p-D-l-thiogalactopyranoside (IPTG). Positive 
clones were tested by PCR using (AC) 10 /(AG) 10 and M13 + / Ml 3" as primers. 
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Sequencing reactions and analysis of 120 positive clones were carried out 
with an ABI 3730x1 DNA sequencer (Applied Biosystems, Foster City, Califor- 
nia, USA), again following the manufacturer's instructions. These sequences 
were analyzed for repeat motif regions of microsatellites using the software 
SSRHunter (Li and Wang, 2005). Of these, 78 clones had microsatellite motifs, 
and primers were designed with Primer3 software (Rozen and Skaletsky, 2000). 
Polymorphisms of all loci with designed primer pairs were assessed with 92 
individuals in two populations of each species from Nangqian in Qinghai Prov- 
ince (population code: NQ) and Hongyuan in Sichuan Province (population 
code: HY), People's Republic of China (Appendix 1). The PCR reactions were 
performed in a 15-jiL reaction volume containing 0.8 uL of template DNA 
(10-100 ng), 1 .5 |iL of lOx buffer, 0. 15 uL of dNTPs (10 mM each), 0.5 jliL of 
each primer (10 mM), 5 U of Taq (TaKaRa Biotechnology Co., Dalian, China), 
and 1 1 .4 uL of ddH 2 0. The PCR cycling profile included an initial step of 5 min 
at 95 °C; followed by 30 cycles of 50 s at 94°C, 50 s at annealing temperature 
for each primer (Table 1), and extension for 30 s at 72°C; followed by a final 
extension step at 72°C for 7 min. PCR products were then electrophoresed by 
QIAxcel Advanced System (QIAGEN, Hilden, Germany). Out of the 60 primer 
pairs, 21 pairs generated amplification products of the expected sizes in which 
16 primer pairs displayed polymorphisms among the populations of the two 
species (Table 1). 

MICRO-CHECKER version 2.2.3 (van Oosterhout et al., 2004) was used 
to assess null alleles and scoring errors. The number of alleles per locus (A), 
observed (H 0 ) and expected heterozygosities (H e ), deviations from Hardy- 
Weinberg equilibrium (HWE), and linkage disequilibrium (LD) between all 
pairs of polymorphic loci were calculated with GENEPOP version 4.0.10 
(Rousset, 2008). Across the two populations of S. alpina, A ranged from three 
to 18, H Q ranged from 0.043 to 0.870, and H e ranged from 0.126 to 0.950. In 
S. mongolica, A ranged from four to 30, H 0 ranged from 0.040 to 1.000, and 
H e ranged from 0.544 to 0.968. Some loci showed significant deviation from 
HWE (Table 2). 



CONCLUSIONS 

The SSR markers developed here are efficient to estimate 
genetic diversity in S. alpina and S. mongolica. Their use at 
larger spatial scales will provide detailed information about 
the distribution of genetic diversity in both species. Fine-scale 
genetic structure studies will enable us to estimate levels of his- 
torical gene flow in these species. Such information is useful for 
building and testing hypotheses on the history of the Qinghai- 
Tibet Plateau in response to climatic and geologic changes. The 
markers are also expected to be helpful in future studies of 
genetic variation and population ecology in these and other 
species in the subfamily Spiraeoideae. 
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Table 1. Characteristics of 16 microsatellite loci developed in Spiraea alpina and S. mongolica. 
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Locus 




Primer sequences (5 '-3') 


Repeat motif 


Fragment size (bp) 


S. alpina 


S. mongolica 


GenBank accession no. 


SA1 


F: 
R: 


ATGGCACGAACTATTGAATG 
GAATGACACGCAATCTATCC 


(GA) 29 


220-292 


52 


53 


KC894821 


SA2 


F: 
R: 


TCCACCACAAGCCCAGTC 
GAGTAACCCAATCACCACAG 


(AG) 17 


128-168 


53 


53 


KC894833 


SA3 


F: 
R: 


GTCAGCGTAATGCGGATG 
CACAAGGCACCACATAGG 


(AG) 20 


234-280 


53 


53 


KC894822 


SA4 


F: 
R: 


GGAGATCGGCTGAAGAAG 
CCACACCCAATCACAACT 


(AG) 16 


116-164 


53 


53 


KC894823 


SA5 


F: 
R: 


TTCACATCCAAGCAGTTCA 
GGTAAGCATCCAAGTCCAT 


(GA) 15 


248-256 


52 


53 


KC894824 


SA6 


F: 
R: 


GGACTTCTGTTCATACCATAC 
GTCAGCGTAATGCGAATG 


(TC) 40 


280-320 


51 


51 


KC894825 


SA7 


F: 
R: 


TCAATCGCACGACAATCC 
TCAAACCTCAAACCCTAAT 


(CT) 8 


120-160 


53 


53 


KC894834 


SA8 


F: 
R: 


GCGTCCAAGACTAATCCA 
ACCGCCTCAGAACTCACC 


(TC) 10 


119-209 


53 


53 


KC894835 


SA9 


F: 
R: 


ACGGACCTGCGGAGAATG 
CGCTCACAAACCCAACTAACA 


(GA) 5 


137-173 


53 


53 


KC894836 


SA10 


F: 
R: 


GTGAAACGAGCTGAAGGT 
GCCATCTTGGAGATATACGA 


(AG) 2V 


331-381 


50 


52 


KC894826 


SA11 


F: 
R: 


CGAGGTGCTTCAAATTACAA 
ATGGAGGAGGATGCTTAGT 


(AG) 30 


214-238 


53 


51 


KC894827 


SA12 


F: 
R: 


CAATC GAAGAGT GAAG AAGAG 
CCTGGTTACTATAGCAATGGA 


(AG) 15 


240-302 


52 


53 


KC894828 


SA13 


F: 
R: 


CTGTTCAACCTACAATCCAA 
CCAGATCCTAACTACTATCGT 


(AG) 30 


172-228 


51 


54 


KC894830 


SA14 


F: 
R: 


AGAATGCTCGCTACCTGC 
GATGTTCGGGCTTGCTAC 


(TG) 8 


180-220 


53 


53 


KC894832 


SA15 


F: 
R: 


GCTTGGACGGATGGAGAT 
TCGCCAGTCTACTTGCTT 


(AG) 18 


160-250 


52 


54 


KC894829 


SA16 


F: 
R: 


CAGAGGGAAGGAGAAGTCA 
TCCGTCCAAAGTCATCGT 


(GA) 20 


186-250 


50 


50 


KC894831 



Note: T a = annealing temperature. 
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Appendix 1 . Locality information for populations of Spiraea alpina and S. mongolica used in the study. The voucher specimens are deposited in the Herbarium of 
the Northwest Institute of Plateau Biology (HNWP), Xining, Qinghai Province, People's Republic of China. 



Species 


Population code 


Population locality 


Voucher no. 


Geographic coordinates 


Altitude (m) 


S. alpina 


NQ 


Nangqian, Qinghai Province, China 


Chensl6037 


31°58'N, 96°30 , E 


4320 




HY 


Hongyuan, Sichuan Province, China 


Chensl6099 


32°46 , N, 102°21'E 


3654 


S. mongolica 


NQ 


Nangqian, Qinghai Province, China 


Chensl6291 


31°58'N, 96°30'E 


4320 




HY 


Hongyuan, Sichuan Province, China 


Chensl6109 


31°58'N, 96°30'E 


3654 



http://www.bioone.org/loi/apps 



4 of 4 



